Localized modes and bistable scattering in nonlinear network junctions 
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We study the properties of junctions created by crossing of N identical branches of linear discrete 
networks. We reveal that for TV > 2 such a junction creates a topological defect and supports 
two types of spatially localized modes. We analyze the wave scattering by the junction defect and 
demonstrate nonzero reflection for any set of parameters. If the junction is nonlinear, it is possible 
to achieve the maximum transmission for any frequency by tuning the intensity of the scattering 
wave. In addition, near the maximum transmission the system shows the bistable behaviour. 



I. INTRODUCTION 

Recently, a variety of different discrete network sys- 
tems have been studied in order to find the optimal ge- 
ometries for enhancing both linear and nonlinear reso- 
nant wave transmission [B El ■ to reveal the conditions 
of the prefect reflection and Fano resonances HSU, as 
well as to analyze the soliton propagation in the discrete 
structures where the role played by the topology of the 
network becomes important ylBlIl- ne of the major 
issues of those studies is to understand an interplay and 
competition of topology and nonlincarity in the dynamics 
and predict new interesting phenomena. 

Many of such waveguide structures can be described 
as discrete linear and nonlinear networks composed of 
straight, bent, and crossed waveguides in which forward 
and backward propagating waves become coupled to each 
other via one or more crossings (or network junctions). 
The well-known systems for realizing these structures are 
photonic-crystal circuits |lCt micro-ring resonator 
structures with more than two coupled channel waveg- 
uides 12], and discrete networks for routing and switch- 
ing of discrete optical solitons @ . The important ques- 
tion is how the waveguide crossing affects the wave prop- 
agation in the network and how we can modify the junc- 
tion transmission by making it nonlinear. 

In this paper, we study the geometry-mediated local- 
ized modes and wave scattering in the structures com- 
posed of N discrete waveguides ('branches') crossing at 
a common point ('junction'), as shown schematically in 
Figs, nfa-d). This system can have direct applications 
to th e phys ics of two-dimensional photonic-crystal cir- 
cuits [lOLIllj . and it describes the transmission properties 
of the V-splitters [see Fig.^b)], the waveguide cross-talk 
effects [see Fig.^c)], and other types of photonic-crystal 
devices. We show that the intersection point of the dis- 
crete waveguides acts as a complex 5-like defect even in 
the absence of any site energy mismatch, i.e. when the 
crossing point is identical to any of the waveguide points. 
Therefore, the waveguide junction can be viewed as a 
topological defect. We show that this new type of de- 
fect gives rise to two stable localized modes (or defect 
modes) at the intersection point, the staggered and un- 
staggered ones, regardless of the amount of nonlinearity 



FIG. 1: Examples of TV-branched discrete network junctions. 



or site energy mismatch at the junction. Next, we study 
the scattering of plane waves by the crossing point, where 
the incoming wave propagates along one of the branches, 
and the waves are scattered to all the branches. In the 
linear regime, we find that for N > 2 the reflection co- 
efficient is always nonvanishing. This fundamental lim- 
itation prevents us from achieving an ideal Y-splitter, 
unless we optimize and engineer the sites near the junc- 
tion [3. For a nonlinear junction of the crossed discrete 
linear waveguides, we reveal the possibility of achieving 
the maximum transmission for almost any frequency of 
the incoming plane wave by tuning the input intensity. 
For some conditions, bistability in the transmission may 
occur near the transmission maximum. To verify our an- 
alytical results obtained for plane waves, we study the 
propagation of a Gaussian pulse across the junction by 
direct numerical simulations. The numerical results agree 
nicely with our analysis. 

The paper is organized as follows. In Sec. II we intro- 
duce our model. Section III is devoted to the study of 
the properties of a linear junction, while Sec. IV considers 
the case of a nonlinear junction and discuss bistability. 
Our numerical results are summarized in Sec. V, while 
Sec. VI concludes the paper. 
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II. MODEL 



We consider the system created by N identical discrete 
linear waveguide (branches) crossed at a common point 
(junction) which we consider to be linear or nonlinear [see 
Figs.^a-d)] . The corresponding waveguide structure can 
be described in a general form by the system of coupled 
discrete nonlinear equations, 



(fe) - A k V 



n > 1 



N 



(1) 



z0o = 53^1 + (e + A|0ol 2 )0o 



fe=i 

where the index k refers to the branch number (fc = 
1,2, ■■-,N). In this model, we endow the junction 
(n = 0) with a linear impurity e and a Kerr nonlinear 
coefficient A. We notice that a similar discrete model 
can be derived for different types of two-dimensional 
photonic-crystal devices based on the Green's function 
approach I n that case, the complex field 4>^ corre- 
sponds to the amplitude of the electric field. 



III. LINEAR JUNCTION 

First, we consider the linear regime when A = 0, and 
all equations Q are linear. We start our analysis by 
investigating the possible localized states in the system 
In order to do that, we look the solutions in the 
well-known form, 



0« = Aexp{- X M^T} , n>0, 



(2) 



where A is the mode amplitude, r is the corresponding 
evolution coordinate (time in a general case, or the lon- 
gitudinal propagation coordinate, for some problems of 
optics), SI is the frequency, and x is the spatial (or tem- 
poral) decay rate of a localized state. After substituting 
this solution into Eq. (Q, we obtain the relations for the 
frequency SI and the decay rate x, 



SI = 2 cosh x 



By denoting y = e 
Eq. as follows 



= (JV-l)e-* + e . (3) 
" x , we rewrite the second equation in 



(N — l)y 2 + eoy -1 = 0, 
which, in general, possesses two solutions 



1 



V 



N - 1 



± 



+ N - 1 



(4) 



(5) 



which are of opposite signs y + > 0, y < 0. As a result, 
there exist two different solutions of our system 



X ± = 



In y + 
\n\y- 



(6) 
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FIG. 2: (Color online) Spectrum of linear waves of the sys- 
tem Q vs. the impurity strength eo for N — 3. In addition 
to the propagating modes u q = 2 cos q, two localized states 
appear, SI . At the critical values these localized states touch 
the band of the propagating modes. Insets show the ampli- 
tude profiles of the corresponding staggered and unstaggered 
modes for eo = 0. 



which describe unstaggered or staggered localized states. 

In the simplest case of a homogeneous system (e = 0) 
there exist always two bound states for N > 2 [see Fig. [2], 



X 



+ = ln(VA - 1) , X ~ = HVN - 1) + i 



(7) 



Therefore, the junction itself constitutes an example of a 
different class of defect: a topological defect. This defect 
supports two localized states whose decay rates x ± are 
inversely proportional to the number of branches N . 

In an inhomogeneous case when eo ^ 0, the decay rate 
of the localized state can change drastically depending 
on the impurity strength. For example, Eq. Q supports 
the solutions y = 1 (or x — 0), when eo = 2 — ./V and 
y = — 1 (or x = * 7r ) J when eo = N — 2. This means that 
in these cases one of the localized states can disappear. 
By looking at the frequency dependence of the localized 
state SI on the impurity strength e see Fig. [2], we see that 
one of the defect modes touches the band of the linear 
spectrum exactly when the corresponding localized state 
disappears. This happens precisely at eo = ±(A — 2). 

Now we analyze the transmission properties of this new 
type of defect. To this end, we consider an incoming 
plane wave propagating in the branch k = 1 and calculate 
the transmitted waves in other branches. We begin by 
imposing the scattering boundary conditions for Eq. Q , 



Ie lqn 

te iqn 



n > 1, k = 1 , 
n > 1, k > 1 , 



(8) 



where uj q = 2 cos q is the the frequency of the incoming 
plane wave. Continuity condition of the plane wave at 
the junction implies I + r = t. On the other hand, the 
equation for </>o can be written in the form, 



U) q t 



IE- iq + re iq + (N- l)te iq + e t 



(9) 
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FIG. 3: (Color online) Linear transmission coefficient for dif- 
ferent values of the impurity strength eo (and N = 3). This 
plot shows that at the critical values eo = ±(2 — N), the 
transmission coefficient becomes maximum at the correspond- 
ing band edge, T(0, 71") = 4/iV 2 , when one of the localized 
states touches the edge of the linear spectrum [see Fig. 
For intermediate values, 2 — N < eo < N — 2, the maximum 
transmission is inside the propagation band. The shaded area 
indicates the possible values of the transmission coefficient. 

By combining these two conditions, we can obtain the 
amplitude of the transmitted wave, 



t = 



2i sin q I 



I 



2ismq + e + {N - 2)e^ JV 1 



(10) 



where a q — (eo + (N — 2) cos q) / (N sin q) . The transmis- 
sion coefficient for any branch is defined as T — \t/I\ 2 . 
From Eq. I|10[l we see that even in the homogeneous sys- 
tem (eo = 0), the waveguide junction produces an effec- 
tive geometry-induced complex (5-like scattering poten- 
tial, e q — (JV — 2)e lq , which depends on the parameters 
of the incoming plane wave. 

We can write the transmission coefficient in the form 



T = 



1 



N 2 l + a 2 n 



(11) 



condition a q — —i, by assuming that the wave number q 
is complex q = q r + i\- Complex wave numbers corre- 
spond to quasi-bound states with the complex frequency 
O, and they can be interpreted as resonances where the 
real part of the frequency gives the resonance frequency 
and its imaginary part describes the lifetime (or width) 
of the resonance. After some algebra, we find that there 
exist only two solutions for q r = 0, 7T, and for \ we have 
exactly the same solutions as in Eq. (J2J. This means 
that our system does not support any quasi-bound states 
(0 < q r < 7r), and there exist only two real bound states 
(q r = 0,n, X ^0). 

What really happens is quite the opposite. Due to the 
frequency dependence of the effective scattering potential 
t q , it simply disappears when the condition eo + Re{e q ) = 
(or a q — 0) is satisfied. And the system becomes al- 
most transparent. The nonzero reflection exists due to 
the nonzero imaginary part of the effective scattering po- 
tential Im(e q ) = (N — 2) sing. 

Equation for poles of the transmission amplitude a q — 
—i differs from the equation for the maximum transmis- 
sion a q = 0, and there is no relation between them. But, 
in fact, there exists the condition when these two equa- 
tions may produce the same results, e.g. when q — 0, n 
and one of the bound state disappears. In that case we 
observe the maximum transmission at the corresponding 
band edge. This result can be understood in terms of 
Levinson theorem, where a bound state just enters to or 
emerges from the propagation band and forms a quasi- 
bound state. 

Now we can compare the dependencies of the defect- 
mode frequencies fi^ = 2 cosh x* and the transmission 
coefficient T on the impurity strength eo [see Fig. 
When eo = 2 — N, the first unstaggered localized state Sl + 
touches the upper band edge, and the maximum trans- 
mission takes place at that edge. For intermediate val- 
ues of the impurity strength, 2 — N < eo < N — 2, we 
observe the maximum transmission inside the propaga- 
tion band which moves towards another edge. Finally, 
when eo = — 2 the second staggered localized state f2 - 
touches the bottom edge together with the occurrence of 
maximum transmission. 



the reflection coefficient is: R = 1 — (JV — 1)T. The 
maximum of the transmission coefficient T max = 4/A 2 
can be reached for a q = 0. We note that i? m i n = 1 — 
4(JV — 1)/JV 2 > for N > 2, and it never vanishes. As 
a result, the JV-junction system described by the model 
(|T|) will always reflect some energy back into the incoming 
branch. This is a fundamental limitation which does not 
allow us to build a perfect F-splitter. 

At this point, the following natural question arises: 
What happens in the system at the maximum transmis- 
sion? First, we notice that this point does not correspond 
to any resonance. According to the linear scattering the- 
ory all (quasi-)bound states of our system can be ob- 
tained as poles of the the transmission amplitude (|10fl . 
The poles can be found in the complex plane from the 



IV. NONLINEAR JUNCTION 

Now we consider the nonlinear junction and analyze 
both the localized states and wave transmission. For 
A ^ 0, the junction can support nonlinear localized states 

. The decay rates x can be found from the equation 
similar to Eq. 



(N - l)y 2 + e~ y -1 = 



(12) 



where we have renormalized the impurity strength, e~o = 
eo + A|j4| 2 . Thus, all previous results about linear lo- 
calized states remain qualitatively valid in the nonlinear 
regime. The only difference is that in the latter case, 
there is an additional dependence on the intensity of the 
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localized state, \A\ 2 . For instance, when nonlinearity is 
focussing A > 0, we have e~ > e . From Eq. (JSJ it fol- 
lows that y + decreases while y~ increases. This, in turn, 
implies that the corresponding modes Q + and f2 _ be- 
come narrower and broader, respectively. In the limit of 
large nonlinearity, the staggered localized mode will ex- 
tend all over the branches, while the unstaggered mode 
will remain confined to essentially one site of the junction. 
Stability analysis shows, that in both cases, the localized 
modes remain stable. We find no other localized modes, 
even in the strong nonlinearity regime. 
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FIG. 4: (Color online) Nonlinear transmission vs. intensity of 
the incoming wave \I\ 2 for different values of the wave number 
q. Other parameters are: TV = 3, eo = — 1, and A = 1. The 
maximum transmission can be achieved for any value of the 
wave number q by tuning the intensity |7| 2 of the incoming 
wave. The shaded area indicates the possible values of the 
transmission coefficient. 

The presence of nonlinear junction generates much 
richer dynamics for the transmission. The transmission 
coefficient can now be written in the form, 

4 1 

T = WJ^+i) ' (13) 

where x satisfies the cubic equation 

(x 2 + l)(x - a q ) - lq = 0, (14) 

j q = A4/ 2 / (TV 3 sing) is the normalized nonlinear pa- 
rameter, and c q = 2 sing. The maximum transmission 
T m ax = 4/7V 2 takes place when x = is a solution of 
Eq. l|14Jl . i.e., when a q = — 7g. It implies that the maxi- 
mum transmission can be achieved for any frequency by 
a proper tuning of the intensity of the incoming wave 
[Fig. [2|. Moreover, the analysis of the cubic equation 
(|H|) reveals that for \a q \ > V3 three solutions are possi- 
ble, and the bistable transmission should occur (FigJSJ. 
We summarize all those scenarios in Fig. [S] We notice 
that the maximum transmission curve, "f q = —a q , lies 
almost at the boundary of the bistability region. 
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FIG. 5: (Color online) Bistable nonlinear transmission in the 
nonlinear junction for q — 2.9. Other parameters are the 
same as in Fig. 2] Bistability takes place near the maximum 
transmission. The shaded area indicates the possible values 
of the transmission coefficient. 

V. RESULTS OF NUMERICAL SIMULATIONS 

In order to verify our theoretical results, we perform di- 
rect numerical simulations of the system under more 
realistic pulse propagation. We launch a Gaussian pulse, 
along the branch k — 1 and study numerically its propa- 
gation through the junction, 

0«(O) = I exp{- ( "~ 2 no)2 - iq (n - n )} , (15) 

where go is th e pulse momentum, / is the maximum am- 
plitude of the wavepacket, a is the spatial width, and hq 
is the initial position. 




FIG. 6: (Color online) Bistability diagram described by the 
solutions of the cubic equation 114H . For the system param- 
eters inside the shaded region there exist three solutions of 
the cubic equation, and the bistable nonlinear transmission 
should be observed (see Fig. 0. The solution x — corre- 
sponds to the maximum transmission T max = 4 /TV 2 , and it is 
shown as a dashed line 7 9 = — a q . 
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(b) |l/ = 6.25 
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FIG. 7: (Color online) Evolution of a Gaussian pulse scattered 
by the N = 3 nonlinear junction (a nonlinear Y-splitter), for 
two different pulse intensities: (a) \I\ 2 — 2.5 x 10 -3 , and (b) 
/ = 6.25. Other parameters are: eo = —1, A = 1, a = 20 and 
go = 2.75. These results should be compared with the results 
for the plane-wave scattering shown in Fig. 0] 

For our numerical simulations, we consider the case of 
a nonlinear V-junction splitter (N = 3) with eo = — 1. 
and A = 1. [see Fig. 0]. Figure shows our results for 
the wave number q$ — 2.75 for two different intensities of 
the Gaussian pulse. In this figure, we show the temporal 
pulse evolution in two branches only (fc = 1 and fc = 2) 
since, because of symmetry, the evolution in the branches 
k = 2 and k = 3 coincide. 

because the evolution in the third branch k — 3 coin- 
cides with that in the branch k — 2 due to a symmetry. 

When the pulse intensity is small [see Fig. 0a)] , the 
whole Gaussian pulse is reflected back into the incom- 
ing branch k = 1, in agreement with our theoretical re- 
sults [Fig.Q]. For larger values of the pulse intensity [see 
Fig. 0b)] , we observe an almost optimal splitting of the 



Gaussian pulse into the branches k — 2 and k = 3 with 
the maximum transmission [see Fig. 0] . We notice that 
in this case the junction 0o remains highly excited even 
after the pulse already passed through it. The excita- 
tion of the nonlinear localized state at the junction is 
possible because of eo = — 1, one of the localized states 
interact strongly with the linear spectrum band (and its 
eigenvalue touches the band edge, see Fig.0). 



VI. CONCLUSIONS 

We have analyzed the linear and nonlinear transmis- 
sion through a junction created by crossing of N identi- 
cal branches of a network of discrete linear waveguides. 
We have revealed that for N > 2 such a junction be- 
haves as an effective topological defect and supports two 
types of spatially localized linear modes. We have stud- 
ied the transmission properties of this junction defect and 
demonstrated analytically that the reflection coefficient 
of the junction never vanishes, i.e. the wave scattering 
is always accompanied by some reflection. We have con- 
sidered the case when the junction defect is nonlinear, 
and studied the nonlinear transmission of such a local 
nonlinear defect. We have demonstrated that nonlinear- 
ity allows achieving the maximum transmission for any 
frequency by tuning the intensity of the incoming wave 
but, in addition, the system can demonstrate bistability 
near the maximum transmission. We have confirmed our 
analytical results by direct numerical simulations. 
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